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Abstract 

The physics of the Mott transition in the anodic region of NiO Resistance Random Access Memory 
(RRAM) is discussed from the Hubbard model. The Hubbard approximation is examined in details and 
it is shown that the Wannier functions in the definition of Hubbard U should not be replaced by atomic 
s-functions when it comes to the metallic side of the transition. The corresponding effective Hubbard U 
is subject to variations, which may also be understood by introducing an effective permittivity of the solid 
as an ansatz. Furthermore, the transition could be demonstrated in the Brinkman-Rice picture. Finally, 
the anodic characteristics of such transition show that it is a local Mott transition. Therefore, the unipolar 
switching NiO RRAM can still have asymmetric I — V curves when a capping layer is inserted, which is 
explained qualitatively by quantum mechanics. 
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The Mott metal-insulator transition (MIT) involves collective localization or de- localization of 
electrons in a solid lli|2j]. Practically, Mott transition is usually achieved by adding pressure or 
doping some ions of different sizes. Nevertheless, none of these two approaches could be dynam- 
ical in an electronic device such that Mott transition can serve as the device physics of a resistive 
switching memory in consumer electronics. On the other hand, there has been an electronic device 
Resistance Random Access Memory (RRAM) that utilizes transition metal oxides (TMO) such as 
NiO [3], which can be subject to a Mott transition theoretically. Rozenberg et al. proposed that the 
boundary domains of a NiO RRAM could suffer from a Mott MIT because of strong correlation 
effects [4]-|(]]. They further pointed out that the RESET transition is due to a lowered occupation 



n 

in the bottom domain [6]. However, the mechanism of such a transition was not discussed in de- 
tails, and they predicted that the insulating region in a unipolar switching NiO RRAM should be 
near the cathode (bottom domain) [6], which is against several experimental findings tA\m- 

In a recent paper [9] we established the analytical current-voltage formulae for carbonyl doped 
NiO RRAM by assuming the densities of states in NiO are different on the metal and insulator 
sides. Nevertheless, the RESET transition is beyond the scope of those quantum transport equa- 
tions. Even if SET seems a phenomenon like tunneling or breakdown, the RESET phenomenon is 
more mysterious since it is an increased applied voltage that causes a transition from ON state to 
OFF state. This is different from negative differential resistance (NDR) phenomenon because no 
orders of resistance change can be observed in NDR. According to this peculiarity, many groups 
have attributed RESET to the thermal rupture of conduction filaments |10]-Jl3], because an in- 
creased voltage leads to stronger Joule heating. Nevertheless, such explanation has difficulty in 
that the power through the device can be higher at the SET point than the RESET point and one 
cannot explain why the filaments are not ruptured right after SET I9Q. Our starting point of the 
RRAM model was the fact that the insulating region in NiO RRAM only lies near the anode, as 
confirmed by the following experiments: (i) the sputtering damage experiment carried out by Ki- 
noshita et al. 020; (ii) the transmission electron spectroscopy (TEM) observation by Park et al. 
[|8|]. These two experiments, especially the first one, strongly confirm that the local region near the 
anode must be emphasized if the RESET mechanism is to be clarified. Note that near the anode, 
the electron quasi-Fermi level is lowered during the transport and thus there should be a dynamical 
electron deficit occurring near the anode. This implies that the Mott criterion Jl^j could play a 
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role, which states that a solid is an insulator if the free electron concentration (nf) is so low that 

n f 1/3 a B < 0.26 . (1) 

Here a B = 4ireh 2 /me 2 is the effective Bohr radius. When it comes to RESET, the applied voltage 
causes strong Fermi level splitting in the two electrodes and the suppression of electron concentra- 
tion near the anode may incur a local Mott transition. This conjecture is attractive in that the Mott 
transition can be found and utilized in real electronic devices. Yet, notwithstanding its successful 
explanation to RESET, the mechanism of such Mott transition deserves further discussions. In the 
present paper, this local Mott transition will be examined from the well-known Hubbard model 
and the Brinkman-Rice picture. 

The Hubbard model is the general model in discussing narrow 3d band TMOs as well as the 
Mott transition. In writing down the Hamiltonian of a solid, it neglects electron-phonon interac- 
tions and the ionic Hamiltonian. Furthermore, the Wannier representation is preferred such that 
the Hamiltonian writes as: 

H = ^ 7 -/'L r ' /" + 9 ^2 V ' a \v a )a< a lv' a ko > ( 2 ) 
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where 



Tij = J dru* (r- ify h(f) u (r - R 3 ) = ^J^ £ (^) 0) 

k 

is the hopping integral, u represents Wannier functions, and f) is the single electron Hamiltonian; 
v (ij; kl) can be expressed in terms of Wannier functions as: 

2 f r u* (fx - MA uj* (f 2 - Rj) oj (r 2 - R t ) u (fi - R k ) 

v (ij; kl) = / / dndf 2 — ^ J - ^ ■ — J - — ^ J - — ^ . (4) 

47re \ri-r 2 \ 

There are many types of integrals in v (ij; kl), which renders the Hamiltonian unsolvable. Hubbard 
in his original work noted that the Wannier functions could be approximately replaced by atomic 
s-functions (</>) when it comes to a very narrow band [15]. Hence, he assumed that only the single- 
center integral v (ii; ii) is significant and should be retained. It has another name as the Hubbard 

U: 

e 2 , , 0* (n - R t ) (j)* (f 2 - Rj) (f 2 - Ri) (n - R k ) 

U = v (ii; ii) = / / df x df 2 — ^ '- ^ ' ^\ '- — ^ . (5) 

4vre yy " |ri-r 2 | 

By neglecting multi-center integrals, one arrives at the Hubbard model. Later, Brinkman and Rice 

ni6fl started from this model and utilized a variational method by Gutzwiller [17], showing that a 



second order MIT could be predicted. Their scenario is that the Hubbard U (force of localization) 
and the tight binding bandwidth W (force of derealization) are competing with each other. A solid 
with half-filled 3d valence band, such as NiO, can be an insulator if the Hubbard U overwhelms the 
bandwidth W. According to the definition of U [Equation (5)], it should be treated as a constant. 
Therefore, the only way to achieve an insulator-to-metal transition in NiO is to add pressure such 
that W would increase. It then implies that a Mott transition is not practical in NiO RRAM. 
Nevertheless, the above discussions are all based on the Hubbard approximation. As pointed 



out by Herring [18], Friedel Ill9h . and Cyrot Il20h . the practical Hubbard U should not be the 



atomic integral, but rather an effective one. We here argue that the Hubbard approximation can 
be divided into two steps: (i) to neglect all the multi-center integrals; (ii) to replace the Wannier 
functions in v (ii; ii) by atomic functions. In writing Equation (5), one applies the whole two-step 
approximation. However, the argument that U is a constant strongly depends on step (ii). If only 
the first step is applied, then 



2 rr u* [n - Ri) u* f 2 - Rj u f 2 - Ri) <p fi - R k 

U = 4— \ \ dndr 2 —± 1 ^ 7 ^ J — ± '- , (6) 

^e JJ \ri-r 2 \ 

whose value is related to the Wannier functions. In case the corresponding energy band is not 



that narrow, e.g. in a metal [27], the Wannier functions may have nothing to do with the atomic 
functions. If NiO is subject to an MIT, then the forms of Wannier functions are supposed to be very 
different on the two sides. Consequently, we prefer one-step Hubbard approximation for MIT. The 
difficulty lies in that the Wannier functions are unknown, otherwise through 

N 



V 3=1 



(7) 

one can derive exactly the Bloch functions (^). In order to be consistent with the Brinkman-Rice 
picture, the following ansatz can be proposed: assume there is an effective permittivity e such 
that 



g2 rp ^^^-Rij^^-Rj^^-R^^n-R^ 

U = - — / / dr\dr2 r-, —, . (8) 
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The permittivity stands for screening effects related to long-range Coulomb interactions, which 
are totally neglected in the original Hubbard model. Such was previously pointed out by Cyrot 
B20[], and Tsuda et al Il22ll . On the other hand, using this ansatz it is natural that U meta i < U insu i at0 r 



because the charge screening is strong when some of the 3d electrons become itinerant. Also, one 



can envisage a first order MIT in NiO. The key parameter is the free electron concentration that 
determines the strength of charge screening. A larger n/ leads to a larger effective permittivity, a 
smaller Hubbard U, and finally a higher conductivity. To sum up, the relations shown in Fig. 1 are 
followed. 



Our argument is similar to that of Mott's 1949 paper [12311 . Mott claimed that once a few 
electrons get delocalized, there is no reason why no more electrons become delocalized because 
the existing itinerant electrons would function. He deduced that the transition should be of first 
order. On the other hand, the present paper claims that no intermediate states exist between the two 
sides of NiO, though NiO could be more metallic due to a further accumulation of electrons. The 
introduced ansatz allows a schematic demonstration of such Mott transition within the Brinkman- 
Rice picture, as in Fig. 2. A traditional way to achieve Mott transition in NiO is by adding pressure, 
which is to increase W and hence the critical value U cr , while keeping U constant. The present 
way, on the contrary, is to reduce U while keeping U cr constant. 

Granted that the metal and insulator sides of NiO could be described from screening effects, 





FIG. 1: The interrelated criteria for a metal and an insulator in a NiO-like TMO, where nj is the free 
electron concentration. 
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another question arises how one can switch from one side to the other. There are two approaches 
to reduce the free electron concentration in NiO: (i) to localize the itinerant 3d electrons; (ii) to 
reduce the total electron concentration. Approach (i) faces a cyclic difficulty that once the 3d 
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FIG. 2: Illustration of the two ways to transform NiO into a metal: (a) by increasing W and thus the critical 
value U cr ; (b) by decreasing the Hubbard U. Here G represents the conductance. 
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FIG. 3: (a) The three regions in NiO RRAM and the average quasi -Fermi level through the device when a 
voltage is applied; (b) demonstration of the ON-state I — V and the RESET phenomenon. 
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electrons get localized one already has a Mott transition, which requires a reduced rif. Approach 
(ii) seems against Kirchhoff's law in the circuits. Nevertheless, in a non-equilibrium problem, 
there can be local electron deficit and excess at the boundaries of the device [shown in Fig. 3(a), 
assuming ballistic transport] [9]. The total electron concentration near the anode is diminished as 
a voltage is applied, since the electron quasi-Fermi level there is lowered. The peculiar behavior of 
RESET, as illustrated in Fig. 3(b), is then explained as a sufficient separation of quasi-Fermi levels 
in the two electrodes. An interesting point here is the cathodic region. Due to the conservation of 
charges, that region must accommodate extra electrons without incurring further changes or phase 
transitions, such that the anodic region could suffer from the same amount of electron deficit. A 
higher free electron concentration tends to make the cathodic region more metallic. As the three 
regions of the device are in series, the resistance change near the cathode has almost no influence 
on the entire device's behavior. 

That the Mott transition is called a local one should be understood within the framework of an 
electronic device. Admittedly, Mott transition is always due to collective electron behaviors rather 
than defective behaviors. Nevertheless, the possibility of such a Mott transition in an electronic 
device strongly depends on the local property of the anodic region. One of its great consequences 
is the asymmetric behavior of the device. The so-called "unipolar" switching RRAM necessarily 
implies a symmetric I — V behavior. However, this symmetry can be broken once a voltage with 
certain polarity has been applied to RESET the device. A good example of this is to insert an insu- 
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FIG. 4: The device structure of a NiO RRAM with a single capping layer, where a\ is the thickness of the 
capping layer. 
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lating capping layer (undoped or less doped NiO) between carbonyl doped NiO and one particular 
electrode |24J], as shown in Fig. 4. In the OFF state, the 1-D potential through such a structure 
depends on the polarity of the last RESET voltage. Figure 5 illustrates the potential profiles of: (a) 
positive polarity and (b) negative polarity. The interesting question is regarding the SET voltage, 
since this is a conventional tunneling problem in quantum mechanics. The transmission coefficient 
T through a 1-D finite potential barrier Vo is [25]: 

-i 



T{£) 
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sinh 2 {2^a) 



(9) 
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FIG. 5: Schematic diagrams of the transport through a NiO RRAM with a capping layer: (a) positive 
polarity; (b) negative polarity. 
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where k is the wave vector, a is the thickness of the barrier, and 

= vSM (10) 

ft 

actually involves the energy of the electron. When the transmission coefficient is high, as is when 
SET occurs, the formula can be simplified as 

In the cases (a) and (b), we may first assume the same energy of the incident electrons. The 
potential barriers of the capping layer and the insulating NiO generated during RESET are assumed 
to be the same: 



1 

1 + k 2 (ai + a 2 ) 



T a = —TT-, ; -s , (12) 



Tb l + k 2 a 2 l + k 2 a 2 ' (13) 

thus 

T a _ 1 | fc 2 aia 2 (fc 2 aia 2 - 2) 

Tb 1 + k 2 (ax + a 2 ) 2 

Since the transmission is supposed to be strong at SET, one must have k 2 a\a 2 < 2, according to 
Equation (11). Consequently, T a < Tb, which implies that one must apply a higher voltage in 
order to SET the device in the case of (a). This has been confirmed experimentally, where the I — 
V curves under opposite polarities were found to be different. The average positive and negative 
Vset values are 2.164V and 1.652V, respectively 

In conclusion, the Mott transition in the anodic region of NiO RRAM has been described 
through the Hubbard model, given that the Hubbard U is treated as an effective one. The screening 
effect renders a first order Mott transition in the solid as claimed by Mott, and a simple ansatz has 
been introduced where the vacuum permittivity in the definition of U is replaced by an effective 
permittivity. The nature of such transition is further shown in the Brinkman-Rice picture. The 
transition is established on a real electronic device structure, since the applied voltage triggers the 
local electron concentration variations inside NiO. In this way, we show that the transition should 
always be a local one. Accordingly, a possible asymmetric device behavior is predicted. The 
proposed mechanism has made the Mott transition practical in real electronic devices. 
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